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Abstract 

The semiclassical collapse of a homogeneous sphere of dust is studied. After iden- 
tifying the independent dynamical variables, the system is canonically quantised and 
coupled equations describing matter (dust) and gravitation are obtained. The condi- 
tions for the validity of the adiabatic (Born-Oppenheimer) and semiclassical approx- 
imations are derived. Further on neglecting back-reaction effects, it is shown that in 
the vicinity of the horizon and inside the dust the Wightman function for a conformal 
scalar field coupled to a monopole emitter is thermal at the characteristic Hawking 
temperature. 
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1 Introduction 



The canonical quantisation of general relativistic isotropic systems carried out in suitably 
chosen variables leads to the whole dynamics being determined by the Hamiltonian con- 
straint [1] of the Arnowitt, Deser and Misner (ADM) construction [2]. Such an approach 
is particularly useful if one wishes to study the semiclassical regime of a system of self 
gravitating matter [3]. In this note we apply the ADM formalism to a model describing 
the collapse of a spherical object leading to a black hole. Further since one of the more 
intriguing aspects of quantum field theory in the curved background generated by a black 
hole is thermal Hawking radiation [4] , we examine how it arises within such an approach. 

In section 2 we briefly review the classical dynamics of an isotropic perfect fluid with 
zero pressure (dust) represented by a scalar field [5, 6]. In particular we observe that 
with a suitable gauge choice for the radial coordinate [6] , the supermomentum constraint 
can be solved at the classical level and one is left with a Hamiltonian whose flow lines are 
the well known inhomogeneous Tolman solutions [7]. These reduce to the more specific 
Oppenheimer-Snyder model [8] if one also assumes homogeneity (subsection 2.1). 

In contrast with a previous treatment [6] we do not impose the temporal gauge (j) = —t 
which amounts to using the (classical) dust degree of freedom as a clock for the system and 
also restricts the solutions to the homogeneous case only. In fact we shall introduce time 
only after having quantised the remaining Hamiltonian constraint obtaining the Wheeler- 
DeWitt equation. This is not a trivial task and in subsection 2.1 we mention briefly the 
main conceptual problems one encounters: the operator ordering, the singularity at the 
origin and the absence of time. 

In order to overcome all of these difficulties, in section 3 we illustrate a decomposition 
scheme for the Wheeler-DeWitt equation of the Oppenheimer-Snyder model in analogy 
with composite systems such as molecules which have two time (or mass) scales. The total 
wave function is factorised into a matter and a gravitational part [3] and on neglecting 
fluctuations (Born-Oppenheimer or adiabatic approximation), matter will follow gravity 
in the semiclassical limit for the latter. 

Finally, in section 4 we introduce a massless scalar field conformally coupled to gravity 
and we allow it to interact locally with a monopole emitter localised near the horizon 
in coincidence with the last shell of the collapsing dust. This is done with the aim of 
obtaining information about the emitted radiation (temperature and intensity) as the last 
shell approaches the horizon. 
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We use units for which c = kBoltzmann = 1, = SttG, the Planck length is then 
= ^fnfi and the Planck mass is rup = Ti/ip. 

2 The isotropic sphere of dust 

The most general three dimensional isotropic metric can be written as: 

e^'^dC^ + e^^dfJ^ , (2.1) 



where dfi^ is the surface element of a unit 2-sphere, v = i^(C,c) and A = A(C, cr) are the 
gravitational degrees of freedom which depend only on the internal time variable a and the 
radial coordinate A real massive scalar field (f) = (t>{C,-,cr) can be related to the dust 4- 
velocity through Ua = (t>,a and Hamilton's equations of motion for the self-gravitating 
sphere of dust can be obtained by variation of the action functional: 



5 = 47r / dadC 



dv d\ dd) 1 
7r„ — + TTA ^ + TT^ ^ - iVT^ - iVi 



(2.2) 



" da ' "'^ da ' da 

with respect to v, A, (/>, their conjugate momenta tTi^, tta and vr^ = —^/\g\fiU^. Here 
/i = a) is the rest mass density of the dust, g is the determinant of the 4-mctric whose 
pull-back on the 3-surface S^. of constant a is given in Eq. (2.1), A'^ and are the lapse 
and radial shift functions. 

In the ADM formalism, N and A^i play the role of Lagrange multipliers for the Hamilto- 
nian density Ti and the radial supermomentum Ti^ respectively, leading to the Hamiltonian 
and radial supermomentum constraints [6, 9]: 
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(2.3) 



= -6-2*^ (V-z/'tt^-A'tta) -e-2'^7r<^</)' = , (2.4) 

where a prime denotes the radial derivative d/d^ and the above are secondary constraints 
in Dirac's language and arise because of the freedom of scaling the time variable a and 
the radial coordinate ( onT,^. 

The residual gauge freedoms allow one to choose both the comoving gauge U°' = S°q, 
that is = and then the scalar field is a function of time only, and the radial gauge 
^/gA-i/ _ ^ which is possible if X'e^"^, considered as a new canonical variable, is a 
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monotonic positive regular function of ( at each value of a. One then has A as the only 
gravitational degree of freedom and on solving the constraints Eqs. (2.3) and (2.4) one 
obtains a final Hamiltonian which does not contain any spatial derivatives so that the 
dynamics for different radii (shells) decouples. This is a consequence of the fact that dust 
is not subject to other interactions except gravity and Birkhoff's theorem ensures that 
each shell at constant C probes only the total mass in its interior [10]. 

After performing a further canonical transformation r = from the Hamiltonian 
equations of motion one obtains a single equation: 

r2 + 2rr + (l-C^) =0 . (2.5) 

Here and in the following a dot will denote a derivative with respect to proper time 
dr = N da. It is also useful to define p = -\/(l — C^)/e, with: 

+1 o<c<i 

C = 1 (2.6) 

-1 C>1- 

In terms of this new radial coordinate the interior of the sphere of dust is given by < p < 
Po, where Pp is the coordinate of the outer surface of the sphere: < Pg < 1 for e = +1, 
> for e = 0, — 1. Of course, for e = the above definition is undefined (it corresponds 
to flat space) and we simply have ({p) = 1 for < p < p^. 

2.1 The homogeneous case and matching with external Schwarzschild 
metric 

The solutions to Eq. (2.5) were obtained long ago [7] and, since the density K//(p, r) = 
F' /{r'^r') is determined only by the arbitrary function F{p), we can restrict our analysis 
to the case in which the dust is homogeneous at each instant of proper time r by taking 
F{p) = 2p^ Kq, where Kq is a constant. The equations of geodesic motion then become: 

r = pKodr,he (rj) = P Kd {rf) 

(2.7) 

r = To ±Kohe{r]) , 



where Tq is another constant and: 



heiv) 



rj — sinry e = +1 

r]^/6 € = (2.8) 

sinh 77 — 77 e = — 1 , 
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with r) a new time variable related to the proper time through drj = ^dr. The density is 
now given by: 




and the line element by: 




One then sees that the parameter rj is the conformal time and the portion of space-time 
filled with dust has the structure of a homogeneous Robertson-Walker manifold [10]. The 
case e = +1 corresponds to a sphere starting from its maximum radius r = 2 p Kq at 77 = tt 
which finally collapses to a point at = while the other two cases (e = 0,-1) start from 
an initial radius at rj = rj^ > and again collapse at r/ = 0. 

We may now consider the exterior metric: the simplest assumption is to have vacuum 
outside the sphere of dust, so that isotropy and Birkhoff's theorem lead to the unique 
Schwarzschild solution [10] which must be matched at the surface r,, of the dust sphere. 
This requires that each shell of dust move along a geodesic in both metrics [10, 11] implying 
that a shell at p < p^^ moves along r = pKciiij) < in a local metric of the Schwarzschild 
type with a mass parameter M^{p) = Kq p"^ / k. The conformal time r]jj{p) at which each 
shell crosses its own Schwarzschild radius is then given by r^{p) = 2M{p) = pKci{rijj), 
so that from the point of view of a distant observer the first horizon to form is the outer 
one (and therefore the only relevant one) located at: 



which is the Schwarzschild radius of the sphere. This whole picture is also known as the 
Oppenheimer-Snyder model [8]. 

One may further express the interior conformal time rj as a function of the external 
Schwarzschild time t for the geodesic trajectory at the surface r = r,,. This is rather 
complicated [10], but we are only interested in the asymptotic relations for r^. In 

such a limit one has: 



r^ = 2M = 2Kop,^ 



(2.11) 



' -2M ln{tan(r7/2) - tan(77^/2)} 



e = +1 



-2M ln(r7 - r]^) 



e = 



(2.12) 



-2M ln{tanh(?7/2) - tanh(r?^/2)} 



e = -1 . 
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Eqs. (2.12) are of course singular at = 77^ since we know that t must diverge when the 
outer shell approaches the horizon, whereas nothing happens to a comoving observer. For 
all the above three cases the equations can be inverted to give: 

r?-r?,:.C,e-*/2^^C, (1-^) , (2.13) 

with C\ a positive constant. The last factor is the square of the redshift for a photon 
escaping from near the black hole. 

We may summarize our results as follows: the Einstein-Hilbert action for the metric 
(2.10) is given by: 

= [dT\Kk^-€K] , (2.14) 



where = Att /o° dpp^ /^l — e is the spatial volume of the sphere and, in order to 
reproduce the solutions Eq. (2.7), the action for the dust must be: 



K J 



dT , (2.15) 

n J 

with a final Hamiltonian constraint: 

Let us note that Eq. (2.16) has been obtained a/ter having imposed the condition that 
matter satisfy the (classical) equation of state for dust (that is a perfect fluid with zero 
pressure) and after having solved for the classical radial momentum constraint Eq. (2.4). 
This in electrodynamics would correspond to a choice of gauge in order to eliminate 
unphysical degrees of freedom. On having identified the independent degrees of freedom 
one may proceed to quantise canonically Eq. (2.16) by replacing t^k ^ — 
However one must also face the following issues: 

1. Operator ordering. The formal mapping does not determine the quantum kinetic 
term uniquely: 

^ ^ irr^^iTT^i , (2.17) 

where a, 6, c, m and n can be arbitrary integers such that a + 6+c = 1 and m+n = 2 
with m and n positive, since in all cases the same classical limit is obtained. 

2. Singularity at the origin. The wave function ^) = 'if{K,(f)) satisfying the quantised 
Wheeler-DeWitt equation must be zero at the origin {K = 0). 
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3. Definition of time. Eq. (2.16) does not contain a time variable, so one must face the 
problem of the introduction of time in order to recover classical motion. 

A consistent scheme solution to the above difficulties which also allows for the quantisation 
of matter is discussed in section 3. 

3 The semiclassical collapse 

In the previous section we briefly reviewed for completeness the classical description of 
the evolution of a sphere of dust. We finally limited ourselves to the homogeneous case 
which we shall further analyse (a generalization to inhomogeneous collapse can be done 
by considering separately shells at different radii [6]). 

As a guide for quantisation one has that in the semiclassical limit one must recover 
classical motion from some qTiantum state of gravity and matter. For our purpose it is 
convenient to consider a system consisting of a portion of a Robertson-Walker space-time 
plus a scalar field (j) = 4>{t) whose total classical action is given by: 



where /x = m/fi is the inverse of the Compton wavelength of the field 4>- On using the 
equations of motion for dust [5] and identifying the density as in Eq. (2.9) the above action 
becomes: 



where K^i is defined in Eq. (2.7). This is in agreement with what has been obtained in 
section 2: Sd = 5^ + S^j where S^^ is given by Eq. (2.14) and Sj^ by Eq. (2.15). 
The Hamiltonian obtained from the action in Eq. (3.1) is: 




(3.1) 




(3.2) 




(3.3) 



where: 



TT 



K 




(3.4) 
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Canonical quantisation (tt^ ~^ t^k — —ifid/dK, iTcj, ^ tc^ = —ihd/d(f)) then leads to the 
Wheeler-DeWitt equation: 



^iK,(j)) =0 . 



(3.5) 



One may now choose an operator ordering in the gravitational kinetic term given by: 

92 1 



(3.6) 



K dK^ K ' 

and following a previously employed procedure [3] (analogous to that used in molecular 
dynamics) we express ^' in the factorized form: 



(3.7) 



which after multiplying on the L.H.S. of Eq. (3.5) by x* ^ind integrating over the matter 
degrees of freedom leads to an equation for the gravitational part: 

K 



1 



(x I x) 

'^^^ {X\dl\x) 7 



(x I x) 



where we have defined a scalar product: 



Further we have set: 



(x I x) 



ijj = exp ■ 



>X*icl>,K)xicl>,K) . 



{dl)i^, (3.8) 



(3.9) 



K 



A{K') dK' \ ^ 



X 



exp 



(3.10) 



A{K')dK'\x, 



with: 



A = - 



.{x\dK\x) 



(x I x) 

If we multiply Eq. (3.8) by x a^id subtract it from Eq. (3.5) we obtain: 



2 / d'ljj \ dx Kfi'^ 



dK dK 



52 



dK^ 



d' 



K 



X 



(3.11) 



(3.12) 



which is the equation for the matter (dust) wave function. Let us emphasize that by a 
suitable choice of operator ordering we have obtained a satisfactory form at the origin 
[K = 0) for the gravitational part {K ip) of the wave function. 
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3.1 The dust wave function 

If we consider the semiclassical (WKB) approximation to the wave function one has: 



dS, 



(3.13) 



dK n dK 

where Seff is the effective action satisfying the Hamilton-Jacobi equation associated with 
the L.H.S. of Eq. (3.8), that is: 

1 



(3.14) 



In such a semiclassical limit ijj will be peaked at K = Kc (classical trajectory). One may 
then define a (conformal) time variable: 



d _ dlnip d 



Further if the R.H.S. of Eq. (3.12) is small: 



d' 



K 



one gets from Eq. (3.12) the Schrodinger equation for a harmonic oscillator: 



dr) ' 



where we have scaled the dynamical phase: 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



and omitted tp while setting K = Kc which is where the semiclassical gravitational wave 
function is peaked and is a solution to the classical equation of motion obtained from 
Eq. (3.14). The solutions are [12]: 



X = X 



( I XiV ) > 



(3.19) 



where Hj^ is the Hermite polynomial of degree N, a normalization factor and: 



(3.20) 



Further: 



2„,2 



Kc{Hm)={N + 1/2) 



Kr 



{N + l/2)hnK^ = {N + l/2)huc . 



(3.21) 
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We may now consider the classical limit for the matter wave function. We introduce: 



Xci(t>,v) 



E 

N 



(2^iV!)V 



00 V4 ^-i (iV+l/2) We (,-,0 ) (^) 



— :j— exp |- 2 7^ [^^^ - <?!>o cos (wc {rj - rjj)f - i 
+7^ ^0 ^^-i sin (wc (r/ - ?7o)) - 7 7^ ^^-o^ sin {2ujc {v - %)) 



(3.22) 



The superposition in Eq. (3.22) corresponds to a minimum wave packet displaced a dis- 
tance (f)o at T] = Tjo and oscillating about (j) = with a frequency tOc which is constant in 
the adiabatic approximation. In the classical limit, ^ ^ 0, one has 7 — > 00 and the wave 
function Eq. (3.22) only has support on the classical trajectory: 



\xciii>,v)f-^s{ct>-Uv)) , 



with: 



Mv) = 4>o cos{uc{ri-%)) . 

The coefficients in Eq. (3.22): 

(7 0o)^ 



-7^ <Poy4: 



N 



(3.23) 



(3.24) 



(3.25) 



(2JV7v!)i/2 

have a maximum at = A^c = 7^ '/'o^/2 and their square obeys a Poisson distribution 
with average value Nc and standard variation Hence, if we just retain the term 

associated with the maximum in the dust wave function {xc — Xnc ) have for the 

expectation value of the matter Hamiltonian for such a state: 



00 /jy ■.2N 



N=0 



(3.26) 



which is just the value of the classical Hamiltonian Kc H^^ for the oscillator in Eq. (3.24): 



1 



(3.27) 



3.2 The gravitational wave function 



The existence of the semiclassical limit implies that the gravitational wave function is 
peaked around some classical trajectory Kc whose Hamilton-Jacobi action is given by 
Seff in Eq. (3.14). Since {H^ ) is independent of Kc the trajectory will be given by one 



9 



of the Kci defined in Eq. (2.7). This implies that Seff must be equal to Sd = + Sj^ 
given in Eqs. (3.2). On comparing Eq. (3.14) with Eq. (3.2) above, one finds that in the 
L.H.S. of Eq. (3.8) one has: 

K,{H^) = Kci{H^) = -KoKa. (3.28) 

K 

Moreover the R.H.S. of Eq. (3.8) must be negligible: 



^Ka{H^). (3.29) 



Eq. (3.28) means that the mean matter Hamiltonian behaves as dust upon neglecting 
fluctuations (Eq. (3.29) and the R.H.S. of Eq. (3.12)) and matter follows gravitation adi- 
abatically (Born-Oppenheimer approximation). 

It is straightforward to derive an expression for ip which satisfies Eq. (3.8) in the 
semiclassical limit and neglecting fluctuations. Suppose that at time rj^ for which dj^Kd = 
(that is r/o = or TT for e = +1 and r/o = for e = 0, —1) the gravitational wave function 
is a gaussian packet centred at the value Kci{r]g) with width b: 

ji(g,.„)=exp{- '^-^y"»»' } . (3.30) 

where for simplicity we have omitted a normalization factor. At a succeeding time 77 ip 
will be given by: 

i>iK,n) = J dK'G.iK -eK„n;K' -eK„rj,)i^{K';rj,) , (3.31) 

where is the Green's function for the (inverted) harmonic oscillator for e = +1(— 1) or 
for a particle moving in a linear potential for e = [13]. One then obtains: 



ip{K,r]) = exp 



X exp |i ^ K, {dlK) (k, ^l-e{dlKY + 2 6 (i^ - e ifj) | 

X expj-^ {K-Kd{ii)f \ , (3.32) 



where he{r]) has been defined in Eq. (2.8) and: 

_b 

k2 (a2^,)2 + 64 (1 - e (a2/l,)2) 



« = - — : : 77TJ2 ■ (3-33) 
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Prom Eq. (3.32) one immediately obtains: 

^ = ' ^^.('y) + = ~^. + o{n') , (3.34) 

in agreement with Eq. (3.13) to leading order in h. We further note that the total grav- 
itational wave function is Kip{K) which is for = as desired: this solution is a 
consequence of our choice of ordering, which however is not relevant in the semiclassical 
limit. 

The absolute value of the gravitational wave function is: 

|^(if,r?)P = exp {-a^ - ifd(^))'} , (3-35) 

If one now considers the limit Ti ^ followed by 6 ^ (the classical point-like limit) the 
classical trajectory is obtained: 

k2|j;|2 

where ( ^ | | ) is the norm of the complete wave function K ip and is given in Ap- 
pendix A. With the above limits interchanged, that is 6 — with h finite, one would 
obtain: 

a^O =^ J^' , ^0 . (3.37) 

However, we expect that when b becomes smaller than the Planck length ip quantum 
gravitational effects (fluctuations) become significant. Hence it is more sensible to consider 
the limit for which one has h ^ ip (corresponding to a minimum size wave packet of the 
order of the Planck length) and then consider Ip 0. For such a case: 

a ~ £p~'^ GO , (3.38) 

which again leads to Eq. (3.36). 

If we also consider the semiclassical limit for matter from Eqs. (3.22) and (3.28) one 
has: 

-KoK^= NcfinKci , (3.39) 

K 

and, taking into account the junction condition Eq. (2.11): 

M^= p^^ NcfiiJL . (3.40) 

We may then conclude that the total mass of the collapsing sphere is proportional to 
the number (large) of dust quanta Nc- The total euclidean action Se (the gravitational 
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entropy) of the system once the black hole has formed is known to be proportional to the 
area of the horizon [14]. On using our result Eq. (3.40) it also follows that: 



(3.41) 



which implies that Se is also quantised in terms of the dust quanta. 

The above apparently differ from the result obtained previously, ~ m^^fn [15], but 
we observe that the appearing in Eq. (3.40) are energy quantum numbers of the dust, 
while the n in Bekenstein's formula comes from quantising the gravitational Hamiltonian. 
In fact, let us consider the case e = +1. Eq. (3.8), using Eqs. (3.28) and (3.29), becomes: 



^. = 0, 



which may be re-written as [16]: 



92 



+ - 



0,n 7 



2k 



where X„ = K — Ko^n- The solutions are: 

V; = i/n(^n/^p)e-^"'/'^- , 



(3.42) 



(3.43) 



(3.44) 



with Hji the Hermite polynomial of degree n and i^o,n^ = ip^ {n + 1/2). It then follows 
that: 



(3.45) 



which is essentially Bekenstein's formula. Let us note that in the other two cases e = 0, — 1 
no analogous discrete quantisation exists. 

3.3 Consistency conditions 

Let us first see how quantisation modifies the geodesic motion since, even if ?/;(K, 77) is 
peaked on the classical value Kd, one has fluctuations around it. Ehrenfest's theorem [17] 
(with the classical limit) will be satisfied if the average position ( f ) will coincide with the 
classical trajectory up to small fluctuations. That is: 



{ ip I I '^) ) 

What one actually obtains is (sec Appendix A for the details): 



(3.46) 



V^a i^ei (3 + 2 a2 Erfc(-a Ka) + 2 (l + KaA e""' 



a 



Vtt (1 + 2 a2 Kci^) Erfc(-a K^) + 2aKa e""' ^-i'' 
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and in the classical point-like limit, Eq. (3.36) or Eq. (3.38), this expression gives exactly 
the classical trajectory: 

p,{k)^p,K,i. (3.48) 

Ehrcnfcst's theorem, as illustrated in Eq. (3.46), implies that fluctuations are negligible. 
In our case this means that one must have: 

A.I<^^>:''^>'l«(/0. (3.49) 

{K) 

For the wave function ■0 one finds (see Appendix A): 



{^\K'^\ ip) 

(3 + 12 a2 Ka^ + 4 K,i^) Erfc(-a Ka) + 2a (s + 2 R^^^ e""' -^^'^ 



2a2 



Vvr (1 + 2 a2 Kd'^'^ Erfc(-a Ka) + 2aKci e""' ^-'^ 

(3.50) 

so that: 



A 1 



a 



2 



, (3.51) 



in the classical point-like limit, Eq. (3.36) or Eq. (3.38). However if one considers the 
limiting procedure used in Eq. (3.37) one has: 

A 37r , , 

. — - 1 ~ 0.18 , (3.52) 



which implies Eq. (3.46) is not satisfied. 

It is also interesting to consider small black holes. Let us then take Kd small at fixed 
h and b ~ £p, one then has: 

37r - 8 1 , , , 

A ~ -—^ -r^ip, 3.53 
4 ^/^^ a 

which of course implies that one must have rg(77) = p^Kdij}) S> ip. Since we are only 
interested in the limit for which Tq approaches r^j from outside, one must have 2M » 



A ~ £p or: 



M^:^mp . (3.54) 



This coincides with what one would expect from naive quantum mechanical considerations. 
Indeed, if X^h = nh/M is the Compton wave length of the black hole, it should be much 
less than its Schwarzschild radius 2 M for the semiclassical approximation to work. This 
corresponds exactly to Eq. (3.54). 
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The remaining consistency conditions are Eqs. (3.29) and (3.16) to which one may 
add the condition for the vaUdity of the adiabatic approximation, that is that the time 
dependence in the matter wave function induced by the slow variable (gravitation) be 
small. This essentially requires that in Eq. (3.12) the second term on the L.H.S. be 
larger than the third (evaluated using Eq. (3.13)). On using the the solutions Eq. (3.19) 
for N = Nc which will approximate classical dust (analogous results are obtained for the 
"vacuum" state iV = for matter) into the R.H.S. of Eq. (3.8) one obtains (see Appendix B 
for more details): 



(iV.^ + AT, + l) < (iVe^ + A^c + l) , (3.55) 

2M/po. The 



where K^i is evaluated at the minimum value we need, that is r^j/ 
condition Eq. (3.29) then becomes: 



2 K,i{H^) 



< 



Analogously, Eq. (3.16) yields: 



2 K,i{H^) 



Nc'' + AT, + 1 

Qnhp^^ Nj' + Nc + l 
64 2A'c + l 
1 , 



Kh Nc' + Nc + l 



(3.56) 



< 



8/iM3 
1 , 



2A^c + 1 
2Nc + l 



(3.57) 



which is essentially the same as Eq. (3.56). 

Finally from Eqs. (3.12) and (3.15) one must have: 

1/2 



n 



AT 2 + AT, + 1 



< 



3 



1 



iV 2 + Nc+l 



2^2 M^M (2A-, + l) 



1 



e = +l 

e = 
e = -1 



< 1 



(3.58) 



which, for the periodic case e = +1, can be understood as the ratio between the periods 
of the gravitational motion and the one for the matter state Xn ! ^^^^ small is 
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a statement of the adiabatic approximation. Since for the semiclassical wave packet in 
Eq. (3.22) Nc ^ I, the above three conditions are essentially the same and, on using 
Eq. (3.40) to express Nc as a function of M, one obtains: 

M » - , (3.59) 
M 

which means that the Schwarzschild radius of the dust must be much larger than the 
Compton wave length of the dust particles, or: 

» ^ . (3.60) 

Since m <C rup for all known elementary particles, the latter expression is a stronger con- 
dition than Eq. (3.54). The fact that the adiabatic approximation is not valid if Eq. (3.59) 
is not satisfied is not surprising. Indeed in such a case fluctuations, corresponding to the 
creation of matter particles, are large (and the evaporation time small [4]). 

4 Coupling with a conformal scalar field 

An observer modelled by a point-like monopole detector localized at a fixed r > and 
coupled to the radiation field can be used to reveal Hawking radiation [4]. Since the 
detector trajectory approaches the uniformly accelerated one in Minkowski space-time as 
r ^ r^, one may take the radiation field in the Unruh vacuum and proceed as in the 
calculation of the Unruh effect [19]. One then obtains that the observer detects a thermal 
flux with a temperature given by: 

T„ = . (4.1) 

For the fiat case (Minkowski space-time) one can also consider a quantum detector de- 
scribed by a gaussian wave packet and one gets the usual Unruh effect for a classical 
point-like limit as in Eq. (3.36) [20]. 

Here we shall examine a different approach. Let us consider the outer dust shell of 
the sphere situated at p = p,, as a quantum system characterized by an internal energy 
spectrum which is related to the number of dust quanta in the shell itself. Prom the 
point of view of a distant observer the position of the shell is not a classical variable but 
is a quantum observable determined by the gravitational wave function describing the 
semiclassical collapse obtained in the previous section. 
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We consider an isotropic massless scalar field ip = <p{p,i]) conformally coupled to 
gravity [18] and to the outer shell of dust. Its Lagrangian density will be given by: 



C - A 



+ fdr [dp f °°dK S{p - p,) Q{p, t) ^{p, r) , (4.2) 

J J Jo {ip \K^\ ip) 

where R is the curvature scalar, Q{p,t) describes a particle (monopole) emitter and the 
factor 6{p — p^^} \KiJj\'^/{ip \ \ ip) forces the interaction to be localized on the outer shell. 
In particular, in the classical point-like limit (Eq. (3.36)) it leads to 6{p — Pq)6{K — Kd) = 
5{r — r-p). The reason we consider an excited monopole (emitter) is to avoid the problem 
of back-reaction on matter and gravity. Further in our case the observer is distant from 
the dust sphere and static, this is the reason we shall use Schwarzschild coordinates for it. 

Let us assume the emitter has a discrete set of internal energy eigenstates described 
by I £^ ) where E is the energy observed by our distant observer. We suppose the emitter 
is initially in an excited state E^^ and decays by emitting quanta of the scalar field (p to 
a state E^ — hu. One may then estimate, using first order perturbation theory, the total 
probability amplitude P{u!,fj) for the emitter to decay in a finite conformal time. It will 
be given by: 

P^u;,f}) = ^ / dr" / dr' e^- (r/",?7') 



2 



h 



2 



rKdirj")drj" r Kd{7j')drj'e''^^'"-''^Dt{ri",v') , (4-3) 



where Q^; = |( -E-o \Q{P(n^)\ EQ — h,ijj) \ is the absolute value of the matrix element of the 
monopole between the states with energy Eq and Eq — Hu!, t" = t{r]") and t' = t{r}') are the 
Schwarzschild times measured by the distant observer expressed in term of the conformal 
times as given in Eq. (2.12), and fj is an upper cut off such that rj^j < fj < 7)^. Further 
during the interval {rj^ , fj) the static monopole emitter near the horizon is immersed in 
the gravitational wave packet associated with the last dust shell and is the Wightman 
function [18] for the isotropic conformal scalar field in the metric Eq. (2.10) evaluated at 
the same spatial point p" = p' = p, but with the Kd in the denominator "smeared" by 
the wave function Kij) of Eq. (3.32): 

\K"^lJ{K",v'T 



Dt{n",r,') ^ / dK" 
Jo 



Jo {i^{ri')\K'^\^{ri')) K"K' ' ^ " ^ 
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where: 



D,{ri",r,') = { 



h 

Stt^ [cos(ry" — Tj' — ie) — 1] 

h 



e = +l 



e = 0,-l . 



(4.5) 



4^2 (■jy// _ ^1 _ I £^2 

In using these propagators one is Umiting oneself to .s-wavcs. However this is known not to 
be too restrictive since there is a potential barrier outside the horizon in the Klein-Gordon 
equation for the scalar field which, in practice, suppresses all higher angular momentum 
modes. On performing the K" , K' integrations one finds: 



KM') Kciiv') Dtiv", V') = R{r}") Riv') A(r?", v') , 



(4.6) 



where (see Appendix A): 

Rirj) = K^iin) 
= 2aKd 



(V>(r?) \K\i:{n)) 

^aKa Erfc(-a Ka) + e""" ^li 



(4.7) 



V7r(l + 2 a2 Kd^) Erfc(-a Kci) + 2a K^i e""" ^li ' 

It is now convenient to change the variable of integration from the conformal time to 
the Schwarzschild time according to Eq. (2.12) and define: 

' At = t"-t' 



T = {t" + t')/2 . 



(4.8) 



Eq. (4.3), with V'(i" = L) = r]'{t' = L) = f) and r}"{t" = 0) = V(*' = 0) = ?7o, then 
becomes: 



P{u,L) 



J, L-2\T-L/2\ 



d{At) e' 



lAt 



Jo "'^ 4M2 

2\T-L/2\-L 

X R{T + At/2) R{T - At/2) D,{At, T) . 



(4.9) 



We now note that is singular for: 

cos(7?" — rj' — ie) = 1 e = +1 

{rj" -r)' -ie)^ = e = 0,-l. 
Using Eq. (2.13) one has near the horizon: 

t" t' T ( At \ 

n" — n' ~ Cg e~ 2M — e~2M = —2 e^^M sinh ( ) , 

V4M/ 



(4.10) 



(4.11) 
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hence on substituting into Eqs. (4.10) in the hmit e — 0, one finds double poles at: 

Atg = iTriqM, (4.12) 

where q is an integer or zero. Apart from these singularities the integrand in Eq. (4.9) is 
an analytic function in the complex variable At. Thus one can close the integration path 
with a semicircle of radius L — 2\T — L/2\ in the upper half complex At plane and apply 
the theorem of residues. In the limit for L — > oo the integral along the semicircle vanishes 
due to the energy exponential and the possible contributions come from the residues alone. 

Let us further examine the above singularities. In terms of the conformal time near 
the horizon one finds: 

(4.13) 

v' — Vh — e~2M . 

It is clear that Eqs. (4.13) have the opposite sign to the one in Eq. (2.13) for q = 2n + 1 
and therefore we must consider only the residues at even q, that is: 



At2n = 8Tr in M . 



(4.14) 



which arc the poles that are associated with a free fall from outside the horizon. The ones 
with odd q correspond to geodesic motion leaving the horizon from inside. Thus for an 
observer outside the horizon the only relevant poles are the ones given in Eq. (4.14). 
The residues to be summed over then are: 

d 



lim 



d{At) 



(At- Ai2n)'e^"^*A(At,T) 



4M^h T_ 

■ CM 



d 



R{T + At/2) R{T - At/2)e 



u)At 



^ d{At) 

where R is given in terms of rj in Eq. (4.7). Further: 



At=At2r, 



d 



d{At) 



R{T + At/2) 



d 



At=At2r, 



d{At) 



R{T - At/2) 



At=At2n 



SO that the expression Eq. (4.15) simplifies to: 



with/3 = (87rM)/n. 

After performing the At integration the probability amplitude becomes: 



P{^^L) ~ / dTR\T)y e-^^^"" . 



(4.15) 



(4.16) 



(4.17) 



(4.18) 
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where N is the largest integer < (L — 2|T — L/2|)/87rM and we have omitted the integral 
along the contour which does not contribute for L large. We observe that in order to 
include at least one pole inside the contour of integration L must be greater than \At2\ = 
Sir M, which implies: 

r/-r7^ ~ < Qe-^'^ , (4.19) 

and if we consider the shell very close to the horizon, that is values of T large in Eqs. (4.18) 
and (4.19), we obtain a probability amplitude per unit time equal to: 

which is a Planck distribution with the usual Hawking temperature = 1/(3 and R{oo) 
is a non zero constant given by Rirj^j). Let us note that the terms omitted in the approxi- 
mations that lead to Eq. (4.20) are mainly concerned with effects at small time T: in fact 
we let L go to infinity and neglected factors of e~^/^^. The corrections one expects when 
relaxing these assumptions are relevant for small values of T and should lead to transient 
effects when the monopole starts radiating or the horizon is forming. 

In the classical point-like limit, Eq. (3.36), — 1 and one recovers the usual field 
theory in the fixed Robertson-Walker-like background. In this last case Eq. (4.20) in the 
large L limit becomes: 

which is equal to Eq. (4.20) apart from the factor of i?^(oo). If we However take the 
alternative limit, Eq. (3.37), that is 6 — >^ 0, one finds that the dust shell decouples from 
the conformal field, since: 

i? ~ a ^ , (4.22) 

and there is no emission: 

P(u;,L)^0. (4.23) 

One may speculate that this effect can be used to eliminate ultra-planckian effects. In 
fact, it is known that, if u is the frequency of the emitted quanta as is measured by a 
distant observer, a fixed observer located near the point of emission at r = will measure 
instead a blue-shifted frequency: 

.> = U.^JLy\, (4.24) 
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and this expression clearly gives u* > nip/h = ip~^ for Tq sufficiently close to r^. In order 
to probe these modes, one must use an emitter (detector) localized in a region smaller 
than oj"^ ~ £p. So one expects that b in the wave function Eq. (3.32) should be less than 
ip for our collapsing shell to couple with conformal quanta of ultra-planckian energies. 
But this would correspond to the limit in Eq. (3.37) which in turn implies Eq. (4.23): 
ultra-planckian energy quanta are not emitted by the monopole. 

A further point worth noting is that R{oo) given by R{r] = rj^) in Eq. (4.7) actually 
depends on e through a (see Eq. (3.33)) but not through Kci{r]^) = 2M/p^^. This of 
course implies that the probability amplitude in Eq. (4.20) depends on the internal dust 
sphere metric. Indeed since: 

h 



o^{Vh) = 

4VPo' (l-ep„2) + 54(i_2ep„2)2 

one would have in the limit b':^ £p, a{'r]^) Kci{ri^) 0: 



1/2 ' 



(4.25) 



vvr 



4M 



1 



4M 1 



1 



e = +l 



/tF bpo 
AM 



(4.26) 



e = -l 



bp, (2p„2 + i) 

that is one has different emission rates depending on the internal geometry of the sphere 
of dust. 

Let us end with a speculation: we considered the matching condition in the classical 
limit, if one wished to consider quantum mechanical corrections one should replace it by 
the continuity of the gravitational wave function inside the dust (Robertson-Walker) with 
that outside (Schwarzschild) . This suggests the results in Eq. (4.26) should not change 
dramatically immediately outside the dust, implying a form of quantum hair leading to 
information on the geometry in the black hole through the intensity of the radiation. 



5 Conclusions 



We have studied the collapse of a self-gravitating homogeneous sphere of dust in the 
context of canonical general relativity by applying the Born-Oppenheimer approach to 
the coupled dust-gravity system. Such an approach allows for the quantization both 
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of matter and gravity and it is only in the semiclassical limit for the latter and in the 
absence of fluctuations that time emerges and one recovers quantum theory for matter 
on a classical background. This has allowed us, in some measure, to estimate quantum 
mechanical effects both for matter and gravity. 

The gravitational wave packets obtained correspond to the quantum mechanical evo- 
lution of a gaussian packet of initial width b (= £p for a minimum wave packet) centred 
on the classical trajectory K^. This has allowed us to examine the conditions for which 
one recovers the usual "point-like" classical collapse which, for example, was necessary to 
implement the matching (or junction) condition. 

Further the fact matter (dust) is also quantised has allowed us to also estimate matter 
fluctuations which lead to a violation of the adiabatic approximation. The validity of the 
approximation not surprisingly requires that the Schwarzschild radius of the dust sphere 
be greater than the Compton wavelength of the matter (dust) particles, indeed should this 
not be the case one would naturally expect fluctuations corresponding to the creation of 
particles by the gravitational field to be important. 

One of the most interesting quantum effects in black hole physics is Hawking radiation 
and in order to to study this we considered a static monopolc, with an internal energy 
spectrum, situated near the horizon which is being formed in coincidence with the col- 
lapse of the last shell. The monopole couples gravitation to a conformal scalar field which 
has support both inside and outside the dust sphere. We found that a thermal Wight- 
man's function corresponding to a Planckian distribution with the Hawking temperature 
appeared in the de-excitation (or excitation) probability of the monopole. 

Again the fact that our last dust shell is described by a wave function allowed us to 
discuss various limits and differing orders of limits. Let us just mention that although 
one always had a thermal Hawking temperature spectrum, the intensity of the radiation 
depended on the internal dust metric. Should the result also be true just outside the dust 
shell (a sort of quantum mechanical "matching" condition) one would have a correction 
to the classical no-hair theorem. 

The most serious drawback in our approach to Hawking radiation is the neglect of 
back-reaction effects. In fact, one would expect the emitted quanta to alter the metric 
outside the dust (since we no longer have vacuum) and also inside, since the outer shell 
loses mass. We hope to return to this and other points. 
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A Useful integrals 



In the text we defined: 



(A.l) 



where ( ip \ | V' ) is the norm of the wave function K-ij;. If one further defines the following 
differential operator: 



D = — 7T h Kci , 

2o? dKci 



then one has: 



Jo 



In particular one finds: 



^\K\i;) = —r V^aKa Erfc(-a K^i) + e"" ^^^^ 
2 L 



(A.2) 



(A.3) 



(A.4) 



{ip\K'^\i>) 



Vtt (l + 2 Kei^) Erfc(-a Kd) + 2aKd e""' (A.5) 



{■4}\K^\i}) 



A a 



^ a Kci (3 + 2 a' Kd'') Erfc(-a Kd) 

+2 {l + a'Kd') e-«'^='' 



(A.6) 



8a5 



Vtt (3 + 12 + Aa^ K^z^) Erfc(-a K^) 
+2aKd (5 + 2a^Kj) e""'-^^ 



(A.7) 



B Evaluation of the Fluctuations 



Prom the definition Eq. (3.10) one has: 



{Xn\qj^2\Xn) {dK))idK (sK^^'^^^ 
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d 



d 



Prom the identity: 



^ y ( XjV I I Xp ) { Xp I Qj^ I XjV , 



y M I Xjv /y ar^ \ Xn I'm ar^^i Xn I ' 



(B.l) 



dK V"^ ' ^^'J dK ' 
on contracting over ( Xl I for L ^ N one obtains: 



(B.2) 



\ Xl \ Qj^ I Xi , 



( Xl I ai^ I Xn ) 

- ^L 



(B.3) 



where one has il^,^ | ) = ^jv I Xiv ) = ^ + 1/2) I Xat )• Further: 



_3 
~2 
3 7i/U 
2 IT 



(B.4) 



where and a are the creation and annihilation operators for the scalar field cf): 

4> = 



2fiK^ 



(B.5) 



On substituting into Eq. (B.l) one finally has: 
^2 , _ ^ 1 



(x^ 



dK^ 



Xn) 



E 



\ XjV I O T/- o TV' I Xn I 



dK ' dK 

9n^ {N + 1){N + 2) g^V^ N{N-1 



(£;,-£;,^,)2 (E^-E,_j2 
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8K^ 



(N^ + N + l) . 



(B.6) 
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